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Abstract: We study nonlinear dynamics of the flux tube between an external 
quark-antiquark pair in J\f = 4 super Yang-Mills theory using the AdS/CFT duality. 
In the gravity side, the flux tube is realized by a fundamental string whose endpoints 
are attached to the AdS boundary. We perturb the endpoints in various ways and 
numerically compute the time evolution of the nonlincarly oscillating string. As a 
result, cusps can form on the string, accompanied by weak turbulence and power 
law behavior in the energy spectrum. When cusps traveling on the string reach the 
boundary, we observe the divergence of the force between the quark and antiquark. 
Minimal amplitude of the perturbation below which cusps do not form is also inves¬ 
tigated. No cusp formation is found when the string moves in all four AdS space 
directions, and in this case an inverse energy cascade follows a direct cascade. 
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1 Introduction 


The gauge/gravity duality ]I|-0 is successfully applied to investigating strongly cou¬ 
pled gauge theories. Through this duality, it is hoped that one can access their 
nontrivial aspects that are hard to be handled because of the strong coupling. In 
particular, among advantages in using gravity duals, it is worth notifying that dy¬ 
namics in time-dependent systems can be powerfully computed from time evolution 
in classical gravity. Applications range over far-from-equilibrium dynamics governed 
by nonlinear equations, and using numerical techniques for solving them attracts 
much attention. For instance, physics of strongly coupled plasma of quarks and glu¬ 
ons at RHIC and LHC brings motivations to numerically study dual gravitational 
dynamics; a series of seminal works is in Refs. [|-|8|. 

Far-from-equilibrium processes in the D3/D7 brane system dual to J\f = 2 super- 
symmetric QCD have been recently studied in Refs. P~|I3[, where partial differential 
equations for time evolution were solved numerically. As a phenomena characteristic 
in non-linear dynamics, it has been found that long time evolution of the D7-brane 
generates a singularity on the brane, and this is understood from the viewpoint of 
weak turbulence on the D7-brane: The energy in the spectrum is transferred from 
large to small scales ITOLO]. Small scale fluctuations there correspond to excited 


heavy mesons in the dual gauge theory. The turbulent behavior of the D7-brane can 
be interpreted as production of many heavy mesons in the dual gauge theory, and 
the singularity formation is interpreted as deconfinement of such mesons. 

To gain a deep insight into this kind of nonlinear dynamics in the gauge/gravity 
duality, in this paper, we will consider a string in AdS dual to the flux tube between 
a quark-antiquark pair in Af = 4 super Yang-Mills theory, and fully solve its nonlin¬ 
ear time evolution with a help of numerical techniques. This setup corresponds to 
focusing on a Yang-Mills flux tube compared with the collective mesons described 
in the D3/D7 system. In addition, working in this setup is simpler than using the 
D3/D7 system and will provide a clear understanding of the turbulent phenomena 
and instabilities in probe branes in the gauge/gravity duality. To initiate the time 
evolution, we will perturb the endpoints of the string for an instant. Our setup is 
schematically illustrated in Fig. [|. The endpoints are forced to move momentarily 
and then brought back to the original locations. We will loosely use the terminology 
“quench” for expressing this process. This action introduces waves propagating on 
the string, and we will be interested in their long time behavior where nonlinear¬ 
ity in the time evolution of the string plays an important role. When we perform 
numerical computations, we will use the method developed in 0, which turned out 
to be efficient for solving the time evolution in probe brane systems. We are also 


motivated by the weakly turbulent instability found by Bizon and Rostworowski [[14 
Our setup may give a simple playground to study that kind of phenomenon. 

The string hanging from the AdS boundary is one of the most typical probes 
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Figure 1. A schematic picture of our setup. Perturbing the endpoints induces fluctuations 
on the string. 


in the gravity dual. This gives the gravity dual description of a Wilson loop corre¬ 
sponding to the potential between a quark and an antiquark |15], |16[. Although the 
potential in a conformal theory is different from that in real QCD, linear confine¬ 
ment can be realized in nonconformal generalization |L7|. In finite temperatures, a 
string extending to the AdS black hole corresponds to a deconhned quark [IB], 0. 
and this has been utilized for studying the behavior of moving quarks in Yang-Mills 
plasma |20|-|26|. Moving quark-antiquark pairs were also considered (27|. [k|. In far- 
from-equilibrium systems, holographic Wilson loops have also been used as probes for 
thermalization [29, BO]. Nevertheless, veiled by these applications to QGP, nonlinear 


(and non-dissipative) dynamics of the probe hanging string in AdS has not been shed 
light on so much. Some analytic solutions of non-linear waves on an extremal surface 
in AdS have been studied in Ref. f31|. Notice that that configuration corresponds to 


a straight string in the Poincare coordinates. 

When it comes to nonlinear dynamics of a string, formation of cusps would be 
primarily thought of. In fact, it is well known in flat space in the context of closed 
cosmic strings that cusp formation is ubiquitous [f32[. Cusp formation of fundamental 
strings ending on D-branes has been also found in Ref. [[h|. We will turn our attention 
to whether there is such formation of cusps also in AdS.[] 

The organization of the rest of this paper is as follows. We start from reviewing 
the static solution and linear perturbations in Section | and where we introduce 
a parametrization convenient for our use. In Section we explain the setup for 
our time dependent computations. We introduce four patterns of quenches that we 
consider, derive the evolution equations and the boundary conditions, prepare initial 


'-A development of a cusp in a decelerating trailing string was discussed in Ref. 
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data, and explain measures for evaluating the time evolution. Sections g,|, and g are 
reserved for numerical results: In Section g, we discuss two of the four quenches where 
the oscillations of the boundary flux tube are restricted to compression waves and 
therefore we call them longitudinal. We evaluate cusp formation, turbulent behavior 
in the energy spectrum, and the forces acting on the quark endpoints. In Section g, 
we show results of a quench where the flux tube oscillates in one of its transverse 
directions. Finally, in Section g, we examine the last quench where the motion of the 
flux tube is in all three spatial directions of the boundary 3+1 dimensions. Section g 
is devoted to summary and discussion. In appendices, we explain details for the 
numerical computations and derive some formulae used in the main text. 


2 A review of the static solution 


We briefly review the holographic calculation of the static quark-antiquark potential 
in the near horizon limit of extremal D3-branes 
AdS 5 x S 5 , 


15|, pT6H ■ The background metric is 


ds 2 * = ~r (—dt 2 + dz 2 + dx 2 ) + £ 2 dQl , ( 2 . 1 ) 

where £ 4 = 4:7ig s N c a' 2 and x = ( xi,X 2 ,x^ ). We consider a rectangular Wilson loop 
with the quark-antiquark separation L along the aq-direction where the quark and 
antiquark are located at X\ = ±L/2. The dynamics of the string is described by 
Nambu-Goto action, 

s = -2 hS dTda ^' (2 - 2) 

where 7 = det( 7 a b) and 7 a b is the induced metric on the string. It is convenient 
to take a static gauge where the worldsheet coordinates (r, a) coincide with target 
space coordinates as (r, a) = ( t,z ). The static solution is then described by a single 
function 27 = Xi(z)Q The Nambu-Goto action becomes 

5 = ~ J dtdz , (2.3) 


where A = 47t g s N c is the’t Hooft coupling. 

Solving the equation of motion of X\{z) gives the bulk string configuration. Let 
z = Zo specify the bulk bottom of the string reached at x — 0, where the regular 
boundary condition d x z = 0 is imposed. The embedding solution is given by 

f 1 w 2 

Xi(z) = ±z 0 dw - = ±z 0 [r 0 + F(z/z 0 -,i) - E(z/z 0 ]i)\ , (2.4) 

Jz/z 0 V 1 - W 

2 We will use small letters for coordinates and capital letters for functions specifying the string 

position. 
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Figure 2. The string profile of the static solution. 


where To = v^r 3 / 2 /T(l/4) 2 ~ 0.599, and F and E are the incomplete elliptic 
integrals of the first and second kinds defined as 


F{x\ k) — I dt 


o v/(l-t 2 )(l -k 2 t 2 ) ’ 


E(x] k) — I dt 


1 — k 2 t 2 
1 -t 2 


(2.5) 


Setting z — 0 in (^.4|), we obtain L/2 = z 0 r 0 . In Fig. ^ we show the prohle of the 
static string in the (z, X\ )-plane. We will use this solution as the initial configuration 
for our time evolution. 

It is also known that the dependence of the potential energy on L is Coulomb. 
The energy is evaluated from the on-shell action, which in general diverges at the 
boundary, but this divergence can be regulated by comparing with the diverging 
energy of two strings straightly extending to the Poincare horizon. The regularized 
energy is then given by 


Eies T(1/4) 4 L • (2 ’ 6) 

In this sense, the quark-antiquark potential in the AdS background does not cor¬ 
respond to the confining potential of real QCD. This, however, is considered as a 
handy playground for testing nonlinear evolution in the gauge/gravity duality. 


3 Linear perturbation theory 

Linearized fluctuations and stability of holographic quark-antiquark potentials have 


been studied in Refs. |35-|38|. Here, we solve the linearized fluctuations of the static 
solution in coordinates convenient for our use. In (|2.4|) , we parametrized the location 
of the string in terms of the ^-coordinate. In this coordinate, however, the static 
solution X\ {z) becomes multi-valued, and this may not be suitable for considering 
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(a) Profiles of special functions (b) Polar-like coordinates 


Figure 3. (a) The profiles of the functions /(0) and g((/>). (b) r = const, and 0 = const, 

surfaces are shown in the (z, xi)-plane. 


linear perturbations. Instead, we introduce new polar-like coordinates (r, 0) in which 
the static embedding is expressed by single-valued functions, 


Z = rf((p) , Xi = rg(<j>) , 


(3.1) 

where we define the functions / and g as 



f (0) = sn(0;i) , 


(3.2) 

M = -{* #'/W 2 = l' # ’* £;(sn( * i);i) + ro 

4/2 \0 + F(sn(0; ? );i)-r o -0o 

(</> < A/2) 
A>A/2) ’ 

(3.3) 


where sn(x; k) is a Jacobi elliptic function defined as the inverse function of F(x] k ) 
given in Eq. (0): F( sn(x; fc); k) = x. For k = i = y/—l, the Jacobi elliptic function 
has roots at x = f3 0 n (n G Z), where (3 0 = 7r/(2r 0 ) ~ 2.622. We find that there is a 
nice relation between / and g, 

Z'O) 2 + g\4>) 2 = 1 . (3.4) 

With these functions / and g , the static embedding (|2.4|) is simply given by r = zq. 
The profiles of the functions /(0) and g(4>) are shown in Fig. |3(a)| . We also depict 
the r = const, and 0 = const, surfaces in the (z, Xi)-plane in Fig. 13(b)) . 

Using t and 0 as the worldsheet coordinates, we can describe the dynamics of 
the string in terms of three functions, 

r = R(t, 0), x 2 = X 2 (t,0), x 3 = X 3 (t,0). (3.5) 

We consider perturbations around the static solution, R = z 0 ,X 2 = X 3 = 0, as 

R(t,4>) = Z 0 {1 + Xi(t,0)} , X 2 (t,4>) = ZqX2 X 3 (t,<j>) = Z 0 X3 (t,4>) , (3.6) 
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Figure 4. Normal mode frequencies of the longitudinal and transverse modes. The 
eigenfunctions for n = 1,2,3,4 are also shown in the inset in each figure. 


where Xi (i — 1,2,3) are dimensionless perturbation variables. We will refer to Xi 
as the longitudinal mode and Xi (* — 2, 3) as the transverse modes. Then, in the 
second order in y, the Nambu-Goto action becomes 


S 



K<t>) ( z o Xi 


xi 2 ) 


E 


its f 2 w 


/ 2 • 2 f 2 \ 

Go Xi - Xi ) 


(3.7) 


where we define ' = d t and ' = d^, and introduce h(4>) = [(g/f)'f] 2 - To derive 
the above expression, we used the relation of / and g and omitted the total 
derivative terms. The equations of motion for (xi,X 2 ,X 3 ) are 

09 1 + u) X i = o, n 

[d 2 + n')xi = o, n' 

Operators % and U! are Hermitian under the inner products 

rP o i 

d0h(0)a(0)/3(0) , (a,py = j d</> , (3.9) 

respectively. We denote the eigenvalues and eigenfunctions of J-L and %' as {cv 2 , e n (</>)} 
and {ui' n 2 , e' n (cj))}, respectively. These are labeled by integers n > 1 in ascending or¬ 
der of u n and u>' n . The eigenfunctions are orthonormalized as (e n , e m ) = 5 nm and 

e > e! V = $ 
c ni u nm- 

It is easy to check the linear stability of the static embedding against the longi¬ 
tudinal perturbation as 

i r^° 

ujl = (e n ,Ue n ) = — dcfi h (d^n) 2 > 0 . (3.10) 

z o Jo 


z oh 


(Jft, fix),-- 




(3.8) 


(i = 2, 3) . 
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In the same way, the stability against the transverse perturbations can be also 
checked, c <j' n 2 > 0. 

We numerically determine normal mode frequencies and eigenfunctions. These 
are plotted in Fig. ff]. The spectra are asymptotically linear: ui n ,u>' n oc n for n —* oo. 
In fact, from the WKB analysis, we can obtain z 0 co n ~ 2T 0 (n + l) for n —* oo [[J3|, - 

Our numerical results are consistent with the WKB approximation. 


4 Non-linear dynamics of fundamental strings 

The main focus of this paper is to study nonlinear dynamics of the string, where we 
make use of numerical techniques for solving the time evolution. In this section, the 
setup for this is prepared. 

4.1 Setup 

We consider AdSs x S 5 Q2.1|) as the background spacetime, and take the static solu¬ 
tion (|2.4f) as the initial configuration. We then consider “quench” on the endpoints of 
the string: We move their positions momentarily and put them back to the original 
positions. A schematic picture of this setup is depicted in Fig. [T]. 

Let us denote the two endpoints of the string as x q (t ) and x q (t), corresponding 
to the locations of the quark and antiquark, respectively. In this paper, we consider 
the following four kinds of quenches on x q {t) and x q (t)\ 


(i) Longitudinal one-sided quench: 

x q (t) = + eLa(t),°,°) , x q (t) = 0, , 

where a(t) is a compactly supported C°° function defined by 

ex P [ 2 (ct - f + 4 )] (0 < * < At ) 


(4.1) 


a(t) = 


0 


(else) 


(4.2) 


The profile of this function is shown in Fig. [5|. The flux tube vibrates in its 
longitudinal direction, and motions are not induced in the transverse directions 
by this quench. Thus, in this case, the motion of the string is restricted in 
(2+l)-dimensions spanned by (t,z,x i). 


(ii) Longitudinal Z 2 -symmetric quench: 


x q {t) = ( - + eLa(t),0,0 


Xq(t) = ( - eLa(t), 0,0 


(4.3) 


We simultaneously quench both endpoints in the opposite directions along the 
flux tube. The string motion induced by this quench is invariant under x\ —> 
—X\ and restricted in the same (2+l)-dimensions as (i). 




Figure 5. A compactly supported C°° function for quench. 


(iii) Transverse linear quench: 


x q {t) 



eLa(t ), 0 



(4.4) 


We shake one of the endpoints along ^-direction. String fluctuations in this 
direction are induced by the quench but not in x 3 -direction. Thus, the string 
oscillates in (3+l)-dimensions spanned by (t,z,xi,x 2 ). 


(iv) Transverse circular quench: 


x q {t) 



, eLa(t), ±cLa/ cc(t)( 1 




(4.5) 


where we choose the upper and lower signs for t < At /2 and t > At/ 2, respec¬ 
tively. The orbit of x q is a circle given by (x 2 — eL/2) 2 + x\ = ( eL/2 ) 2 . String 
fluctuations along both x 2 - and x 3 -directions are induced by this quench, and 
hence the string moves in all (4+l)-dimensions spanned by (f, z, X\, x 2 , x 3 ). 


In Fig. ||, we show schematic pictures of the quenches (i)-(iv). These patterns are 
chosen to represent typical string motions, particularly with different dimensionality. 
String dynamics is specified by two parameters e and At once we choose a quench 
type. 

In this paper, we will take modest values for the amplitude of the quench (e ~ 
0.01) since we are interested in nonlinear evolution starting from small deviation 
from the linear theory and focus on weak turbulence and cusp formation driven by 
the nonlinearity. For a large value of e, we expect some effect similar to the overeager 
effect found in Ref. |J: The string will be able to plunge into the Poincare horizon 
z = oo because of the strong perturbation. Such strong quenches will be studied in 
detail elsewhere |39| . 
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Longitudinal Longitudinal Transverse 

one-sided Z 2 - symmetric linear quench 

quench quench 



Transverse 
circular quench 


Figure 6. Schematic pictures of the quenches we will consider. 


4.2 Basic equations 

To calculate the time evolution on the string worldsheet, we find it convenient to use 
double null coordinates. With worldsheet coordinates (u,v), the string position is 
parametrized as 


t = T(u,v ) , £ = Z(u,v ) , x = X(u,v) . (4.6) 

Using these expressions with Eq. 0 . we obtain the induced metric as 


7 uu — t^(—T 2 u + Z 2 U + X 2 U ) , 7 vv — t^{—T 2 v + Z 2 V + X 2 V ) , 

e 2 

luv — ~^{~T u T jV + Z U Z, V + X u ■ X iV ) . 


(4.7) 


The reparametrization freedom of the worldsheet coordinates allows us to impose 
the double null condition on the induced metric as 


Cl — 7 uu 0 , ( '2 — 7mi 0 . 


(4.8) 


Notice that these conditions do not fix the coordinates completely: There are residual 
coordinate freedoms, 

u = u(u) , v = v(v) . (4.9) 

These will be fixed by boundary conditions and initial data. 

In the double null coordinates, the Nambn-Goto action (^.2|) becomes 


S = 


1 



1 


I'uul'vv 


27rcr' 


dudv^f u 


27rcr' 

? 2 r 1 

dudv -^(—T’uT'V + Z U Z V + X u ■ X 


(4.10) 


2-Ka’ 


Z 2 


J ,V 


U 5 
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where at the second equality we eliminate the square root in the action using the 
double null conditions (|4.8| ). (Note that j uv is negative.) From the action, we obtain 
the evolution equations of the string as 


T }UV — -^{T u Z tV + T v Z tU ) , 

Z, uv = —{T, u T jV + — X jU ■ X iV ) , 


-X u „ — — (X iU Z )V + X jV Z tU ) 


(4.11) 


Using these, we find that the constraints (|4.8|) are preserved under the time evolution: 

d v C\ = d u C 2 = 0 . (4.12) 


Hence if we impose C\ = C 2 = 0 on the initial surface and the boundaries, the 
constraints ( |4.8| ) are automatically satisfied in the whole computational domain. 

Our numerical method for solving the evolution equations is briefly summarized 
in Appendix |A] and its numerical error is estimated in Appendix [B]. For more details 
of the numerical method, see also Appendix A in [0. 

4.3 Boundary conditions 

On the worldsheet, there are two time-like boundaries that correspond to the two 
endpoints of the string attaching on the AdS boundary. We need boundary conditions 
there. Using the residual coordinate freedoms ( |4.9|) , we can fix the locations of the 
boundaries to u = v and u = v + /3 0 . The boundary conditions for the spatial parts 
of the target space coordinates are given by 


Z\ u =v 0 , -X|u=i> Xq , Z\ u=v _^_fj 0 0 , X\ u=v _i r p 0 Xq . 


(4.13) 


We also need boundary conditions for T. To derive them, we solve Eq. ([4.11|) and 
(|4~8l) near the boundaries. Defining r = u + v and a = u — v (r e (— 00 , 00 ), 
a G [0,/3q]), we obtain the asymptotic solutions around a = 0 as 


T = to(r) + 

ry t0 . 

Z = —a + 

7 


X = x q (r) 


~io ~ l 2 v ■ ai 2 0 


cr 2 + v ■ x 3 a 3 + 


t 0 _ 1 
67 ~ 2 
~1 


yu • ai 0 t 0 - J( 37 2 a 2 + v ■ j)t 3 0 a 3 + 
o 


(4.14) 


-vt 0 - 7 - 


at 


12 


a 2 + x 3 (r)a 3 H- 


where' = d/dr, v = x q /t 0 , a = v/t 0 , j = a/t 0 and 7 = l/\/l — v 2 . The asymptotic 
solutions near a = (3o can be given by replacing a —> — a and x q — > x q in the above 
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expressions. From the second equation in Eq. (|4.14|) , we have d a Z\ (7=0 = 
Using this, we obtain 



Tr\*= 0 = sJ(Z, a \ a=o y + xl , T T \ a= p 0 = yJ(Z I(T \a=p 0 ) 2 + . (4.15) 

These equations determine the time evolution of T at the boundaries a = 0, /?o- 
Their numerical implementation is explained in Appendix [A], 

For consistency, the speed of the quark endpoints during the quench should be 
slower than light, |u| < 1. Otherwise, the Lorentz factor 7 becomes imaginary. 
Solving this condition, we obtain constraints for the quenches (i)-(iii), 


eL 



(19 - 5\/l3)e v/ ^ _3 
12vWl3-21 


~ 0.1865 , 


(4.16) 


and for the transverse circular quench (iv), 

eL ^ 1 

At < 4^2 ' 


0.1768 . 


The parameter values examined in this paper satisfy these conditions. 


(4.17) 


4.4 Initial data 


Before the quenches are applied, we assume that the string is static, namely, we use 
the static solution ( |2.4|) as the initial configuration. For numerical computations, 
we need initial data written in the (u, u)-coordinates. Substituting Eq. 0 into 
Eqs. ( |4.8|) and ( |4.11|) , we can express the static solution in terms of (u,v) as 

T(u,v) = Zo[<j>i(u) + M V )\ , Z(u,v) =z 0 f((/) 1 (u) ~M V )) i ^ 41g ^ 

Xi(u,v) = z 0 g((t>i(u) - 4>2 {v)) , X 2 (u,v) = X 3 (u,v) = 0 , 


where (j) 1 and (j) 2 are arbitrary functions associated with the residual coordinate 
freedom ( |4.9|) , and the functions / and g are defined in Eqs. (|3.2|) and (|i-3|) . 

Locating the initial surface at v = 0 on the worldsheet, we parametrize the initial 
configuration as 


T(u, 0) = z 0 u , Z(u, 0) = f{u) , 1Q 

Xi(u, 0) = z 0 g(u) , X 2 (u, 0) = X 3 (u, 0) = 0 , 

where we set the free functions 4>i(u) = u and (f) 2 ( 0 ) = 0 so that the boundaries are 
at u — 0, f3(). If the string endpoints are not perturbed, our numerical calculations 
describe the static evolution of the exact solution ( |4.18|) . 


4.5 Quantities for evaluation 

In the non-linear dynamics of the string, we expect to observe formation of cusps. 
For detailed analyses, we will in particular use the following quantities for evaluation. 
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4.5.1 Cusp formation 

Let us consider the profile of the string on a surface with t ^constant. There is 
a cusp if the string does not change its target space position when the parameter 
on the string is varied. Using the w-coordinate as a parameter on the string in the 
t ^constant surface, we obtain the conditions for the cusp formation as d u Xj |i =const = 
0 where Xj = ( Z,Xi ) and i — 1, 2, 3. These conditions are rewritten as 

Ji = T, u Xj )V — T v Xj >u = 0 . (4.20) 

In onr numerical computations, we monitor the roots of J/, which typically appear as 
curves on the ( u , u)-plane. If these curves overlap at a point, we find cusp formation, 
and if the overlap continues in the time evolution, it is implied that the cusps continue 
to exist. 

As an obvious corollary, if (f4.20| ) is satisfied, there is a necessary condition that 

Xi, u Xj tV — X I)V X J)U = 0 (4.21) 


is also satisfied. This condition is conveniently utilized for a consistency check of the 
cusp formation detected by (|4.20|) . 


4.5.2 Energy spectrum in the non-linear theory 

We will also study the energy spectrum of the non-linear fluctuations of the string be¬ 
cause, from the spectrum, it is expected to find weak turbulence on the string. Once 
a dynamical solution (T(w, v ), Z(u, v ), X(u, v )) is calculated, we can convert it to the 
polar-like coordinates introduced in Eq. (0: r = R(u,v ) and 0 = $(«, u).0 Elim¬ 
inating the worldsheet coordinates (u,v) from ( T(u,v ), R(u,v), &(u,v), X 2 (u,v), 
X 3 (u, v )), we can express the dynamical solution using target space coordinates ( t , 0) 
as 

r = R(t, 0), x 2 = X 2 (t, 0), x 3 = X 3 (t, 0) . (4.22) 


As in Eq. (|3.6|) , we define the non-linear version of the “perturbation” variables 
Xu X 2 , X3 as 


XT (0 0 ) =- 1 


z o 


A ^ Xi(t, 0 ) ^ ^ 

Xi\t,4>) = —r- (* = 2,3) 


zo 


(4.23) 


We then decompose yq, \ 2 and y 3 with the eigenfunctions of the linear theory e n (0) 
and e(j(0), which were introduced below Eq. Q3.81) , as 


oo oo 

Xi = Y1 c nOK(0) > Xi = Yl <(^^(0) (*' = 2, 3) . (4.24) 

n =1 n= 1 

3 Solving Z(u,v)/Xi(u,v) = f(4>)/g((l)) for 0, we have 0 = $(u,i;). Then, we can obtain R(u,v) 
from R(u,v) = Z(u, v)). 
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Using the mode coefficients c n and c l n , we define the energy contribution from the 
n-th mode e n (t) and the total energy in terms of the linear theory e(t) as 


En {/'} 


VX. 


zo 


47r 


+ u n c n 


E< 

i=2,3 


xi 2 I . f 2 i 2\ 

I - ^n. ^-'n ) 


e{t) = s ^2,E n . (4.25) 

n =1 


The quantify e n is conserved in linear theory. Therefore, if we find time dependence 
in e n , it is a fully non-linear effect. Note that the total energy £ defined in the linear 
theory is also time dependent in the non-linear theory although the time dependence 
is suppressed by the amplitude of the quench: e/e = (9(e). Since we consider only 
small e in this paper (e ~ 0 . 01 ), we do not put emphasis on the time dependence of 
the total energy. In our actual numerical calculations, we take a cutoff at n = 50 for 
evaluating e in Eq. (|4.25|) . Its cutoff dependence is also not essential for our following 
arguments on the energy spectrum on each f-slice. 


4.5.3 Forces acting on the heavy quarks 

From the dynamical solution of the string, we can read off the time dependence of 
the forces acting on the quark and antiquark located at x = x q and x q . The force 
acting on the quark can be evaluated from the on-shell Nambu-Goto action as 


( F i(t)) = ^7 1 (Sij + 'y Z v i v j )d%X j \ g=0 , 


(4.26) 


where 7 = 1/ yjl — ( dx q /dt ) 2 . The same formula can be applied to the force acting on 
the antiquark {F(t)) by replacing x q —y x q . The derivation of this expression is sum¬ 
marized in Appendix |D]. For our numerical analysis, it is convenient to rewrite this 
expression in terms of the (r, u)-coordinates. Using the asymptotic expansions ( 4.14 ), 
we obtain 


(F(t)) = 


3v/A 


7 2jc 3 


2tt 


i 3 0 


(4.27) 


We need to extract the third order coefficient x 3 from our numerical data. For this 
purpose, it is convenient to define 


Y = X 


X n 


(T - t 0 )v - 7 2 


7 [v ■ a)v 


7 


\ ]a 


(4.28) 


Using Eq. (|4.14| ), we find that this function behaves as Y ~ [*3 — (v • x^v]^ near 
the boundary. To evaluate x%, we fit numerical data for Y by a function ccr 3 in 
a G [0,0.1], 


5 Results for the longitudinal quenches 

From this section to Section [/]. we discuss results of numerical computations for our 
quenches (i)-(iv). In this section, we firstly treat the longitudinal one-sided quench 
(i) and then discuss the longitudinal Z 2 -symmetric quench (ii). 
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5.1 Cusp formation 

We start from the longitudinal one-sided quench (Cl). In Fig. [/], we show snapshots 
of the time evolution of the string for A t/L = 2 and e = 0.03. The left panel is 
just after the quench, t/L = 2, 2.2, 2.4. We observe that perturbations are induced 
on the string by the quench. Late time behavior is shown in the right panel for 
t/L = 7, 7.2, 7.4, where it is seen that cusps are formed on the string. The time for 
the cusp formation is evaluated as t/L ~ 5 by using (|4.20| ), while the snapshots are 
taken when it is easy to confirm the cusps visually. Although the solution plotted in 
the target space (z, aq)-plane looks singular on top of the cusps, the fields T(u,v), 
Z(u,v) and X\{u,v) are indeed smooth on the (u,n)-plane. Therefore, even after 
the formation of the cusps, the time evolution can continue without a breakdown of 
numerical computations. Physically, however, finit e-N c effects can be important at 
cusp singularities, and the time evolution after the cusp formation may be considered 
as unphysical without such corrections. We will discuss possible finite- N c effects at 
cusps in section ||. For the present, the dynamics of cusps are discussed without 
taking into account these corrections. We find that the cusps appear as a pair. (The 
cusps at t/L = 7.0 are magnified in the inset of the right panel.) The propagating 
speeds of the two cusps are different, and they continue to propagate on the string. 

Following this example, we survey cusp formation by changing quench parame¬ 
ters using the method described in Section |4.5.1| . Results are shown in Fig. || when 
A t/L = 2 is fixed and e is varied. In the left panel, the times for the cusp formation 
for each e are plotted, and the corresponding locations in the ^-coordinate are shown 
in the right panel. The time for the cusp formation becomes longer as the quench 
amplitude becomes smaller. It is also seen that the cusp formation times are mildly 
discretized, as well as the locations of the formation away from the boundary. This 
discretization indicates that the cusps might be difficult to form near the boundary. 

The times and locations of the cusp formation tend to degenerate in the very 
early time before the first reflection of the initial perturbation wave at X = —L/2 
around t/L ~ 3.5. Features in this region seem to be influenced by the largeness of 
the quench and would be different from late time dynamics. Quenches with large 
amplitudes will be reported elsewhere [[39] . 


The cusp formation time becomes longer as the amplitude gets smaller. A natural 
question then is if the cusps can be formed for any small amplitude e. As e decreases, 
however, the variations in the fields are tinier and tinier, and eventually it becomes 
numerically difficult to use the cusp condition ( |4.20|) for finding the cusp formation. 
To obtain a reasonable inference, we fit results and extrapolate to smaller e. A result 
is given in Fig. |D], where we fit data points^ at e > 0.01 and t cusp /L > 10 by a 


4 The points marked with the purple boxes were not derived from directly computing Eq. ( |4.20| ) 
because of difficulty in small e. We are, however, able to find a suspect of cusp formation by looking 
at the plot of the necessary condition ( 4.21 ). We do not use these points for the fit, but they look 
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Figure 7. Snapshots of the string time evolution for the longitudinal one-sided quench 
with A t/L = 2, 6 = 0.03. 



(a) Cusp formation time 
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(b) Corresponding formation points 


Figure 8. Cusp formation for various e when At/L = 2. The formation points in (b) 
correspond to those at the same time in (a). In the plots, the points are computed at every 
0.001 variation of e for t cusp /L < 8 and every 0.0005 for t cusp /L > 8. 


polynomial^ ae + be 1 / 2 + c. From the extrapolation, we find that there is a critical 
value e crit below which the cusp formation time would be infinity. In Fig. we obtain 
Ccrit ^ 7.5 x 10" 3 . 

We repeat this procedure to estimate the critical e for different A t/L and see 
how it changes. For each At/L, we fit data by ae + fee 1 / 2 + c and extrapolate to 
the limit of infinite cusp formation time to read off the value of critical e. Results 
are shown in Fig. 0. We find that the critical value scales as (A t/L) 3 . A fit of our 
results is e cr i t = 9.3 x 10” 4 (Afe/L) 3 . Note that this scaling may be altered if At/L 
becomes very long and the wavelength of the induced wave is comparable with the 

consistent with the fit. 

5 This fitting function is chosen by our intuition, and other choices for extrapolation could be 
utilized. For instance, the data can be also fit with ae 2 + be + c, and a qualitatively consistent 
result can be obtained. 
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Figure 9. Cusp formation in small e for A t/L = 2. The dashed line is a fit of the red 
points. The points marked with the boxes are not used for the fit but are consistent with 
the fit curve. The fit curve reaches (t cusp /L)" 1 = 0 at e ~ 7.5 x 10~ 3 . 
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Figure 10. The dynamical phase diagram for cusp formation. The red points are 
estimated values of critical e at different A t/L. The dashed curve is a fit given by 
e = 9.3 x 10~ 4 (At/L) 3 . Cusps are formed above this curve, but not below. 
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length of the hanging string. The critical amplitude in such a case, however, will be 
also big. For this reason, we do not focus on larger A t/L in this paper. 

The cusp formation found here is similar to the weakly turbulent instability in 
the global AdS [|14|: Small fluctuations in that AdS propagate between the boundary 
and the center and, eventually, the perturbations collapse into a black hole after 
several bounces. A difference between our cusp formation and that AdS instability 
is their critical amplitude of initial perturbations. The AdS instability occurs for 
arbitrary small initial perturbations while our cusps are formed only for e > e cr i t > 0. 
By a perturbative study in Ref. |h|, it has been suggested that a commensurable 
spectrum is a necessary condition for the collapse by arbitrary small perturbations. 
As in Fig. [f], the linear spectrum of the string, u n and u}' n , are commensurable only 
for n —y oo. Because of the non-commensurable spectrum of the string, we need 
finite perturbations for the cusp formation. 
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5.2 Energy spectrum in the non-linear theory 


Given the formation of cusps, we look into the time dependence in the energy spec¬ 
trum following the procedure described in Section |4.5.2| . For samples, we focus on the 
following four cases: (a) At/L = 2, e = 0.005, (b) At/L = 2, e = 0.01, (c) At/L = 2, 
e = 0.03, and (d) At/L = 4, e = 0.07. With the parameters (a), cusps do not form on 
the string, while cusps are created for (b), (c) and (d). In Fig. [11], we show the time 
dependence of the energy spectra for these parameters. The dashed curves are the 
energy spectra computed in the linear theory; see Appendix [C] for the calculations. 
Although the spectra are defined only for integer n, for visibility of the plot these 
results are generalized to continuous n. 

High frequencies in the energy spectra are suppressed when the cusps are not 
formed. In Fig. |11 (a)|, the spectrum can be well approximated by the linear theory 
just after the quench {t/L = 2). Although it slightly deviates from the linear theory 
as the time increases, we do not find any remarkable change in the late time. 

In contrast, the energy spectra show power law behaviors in the cases of cusp for¬ 
mation. In Fig. |ll(b)| , although the spectrum can be well approximated by the linear 
theory just after the quench, we see the growth of the spectrum in high frequencies as 
time passes, apparently because of the nonlinearity in the time evolution equations. 
In Fig. 11(c) and 11(d) , the spectra deviate from those in the linear theory even at 
t = At, and this indicates that the nonlinearity evolves even in the quenching time, 
0 < t < At. For these three cases, we find a direct energy cascade: The energy 
is transferred to higher n-modes during the time evolution. Eventually, power law 
spectra are observed, and these behaviors persist until the time of cusp formation. 
(See magenta points. )[] In particular, as seen in |ll(b)|, the time for reaching the 
power law behavior can be rather earlier than that for the cusp formation, and once 
realized, the behavior lasts until the cusps are formed.[] Hence, the cusp formation 
on the string can be regarded as a variation of weak turbulence [|T1[] ,|j 

We fit the power law spectra by e n oc n~ a . Just before the cusp formation, we 
obtain a = 1.313 ± 0.171, a = 1.354 ± 0.128 and a = 1.432 ± 0.017 for (b), (c) and 
(d), respectively. The exponents are distributed around a ~ 1.4. 


5.3 Time-dependence of the forces acting on the heavy quarks 


Following the procedure in Section [4.5.3| , we compute the forces acting on the quark 
and antiquark as functions of time. Since the string motion is now restricted in the 


6 After the cusp formation, the function R(t, (f>) becomes multi-valued and the energy spectrum 
is ill-defined. Therefore, we only show spectra before the cusp formation. 

' The string is smooth before the cusp formation, and the energy spectrum e„ must fall off faster 
than any power law function as n —> oo. Hence, it is indicated that the power law spectrum will 
be no longer maintained at n> 50, while many numerical efforts are necessary for computing the 
spectrum in such a region. 

s The AdS weak turbulence was found to be characterized by a Kolmogorov-Zakharovscaling ]40|. 
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(a) A t/L = 2, e = 0.005 (no cusp) 


(b) A t/L = 2, e = 0.01 (cusp) 


i 

o.i 


00 

g 

10 



(c) At/L = 2, e = 0.03 (cusp) 


(d) At/L = 4, e = 0.07 (cusp) 


Figure 11. Time dependence of the energy spectrum. For the parameters (a), cusps 
do not form on the string. For the parameters (b), (c) and (d), cusps are created on the 
string. The dashed curves are the energy spectra computed in the linear theory. Red 
points are spectra just after the quench. Magenta points correspond to the time for cusp 
formation. For (b), (c) and (d), the energy is transferred to higher modes as time increases, 
and eventually power law spectra are reached. In (b), light blue points indicate the time 
when the spectrum realizes the power law, which lasts until the cusp formation. We fit the 
magenta points by e oc n~ a , and the results are plotted with dotted lines. 


(z, Xi)-plane, only the aq-components of the forces are non-zero. In Fig. [L2|, we plot 
the forces (Fi(t)) and (Fi(t)) for e = 0.005 and A t/L = 2. For these parameters, 
cusps do not form on the string. Figure 12 (a)| is for the early stage in the time 
evolution (0 < t/L < 20), and Fig. |12(b)| for a long period (0 < t/L < 160). In 
Fig. |12(a)| , pulse-like oscillations are repeated at regular time intervals. The period 
corresponds to the timescale for the fluctuation induced by the quench ( |4.1|) to go 
back and forth between the two boundaries. Because of the dispersive spectrum 
in the linear perturbation, the initially localized oscillations tend to spread as time 
increases as in Fig. |12(b)|, and there is no way for them to converge again. We can 
also see that (F\ (t)) < 0 and (F\ (t)) > 0 throughout the time evolution, and this 
implies that the force between the quark and antiquark is always attractive. 
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(a) 0 <t/L< 20 


(b) 0 < t/L < 160 


Figure 12. Time dependence of the forces acting on the quark and the antiquark for the 
longitudinal one-sided quench with e = 5.0 x 10 -3 and A t/L = 2. In this case, cusps do 
not form on the string. 


For the parameters where the cusps form on the string, the time evolution of the 
forces is different from the previous example. In Fig. [13], we show (Fi(t)) and (Fi(t)) 
for e = 1.0 x 10 -2 and A t/L = 2, where figures (a) and (b) are for 0 < t/L < 13 
and for 0 < t/L < 30, respectively. In figure (b), we take the absolute values of 
the forces, and plot in the log scale in the vertical axis. We find that the pulse-like 
oscillations are getting sharp and amplified as the time increases. The forces can 
change the sign because of their large oscillations, and this implies that the force 
between the quark and antiquark can be repulsive temporarily. Eventually, when a 
cusp arrives at the boundary after the cusp formation, the force diverges. 

A mathematical explanation why the forces diverge after the cusp formation 
is given as follows. Near the AdS boundary, the time dependent solution is well 
approximated by Eq. ( 4.14Q . The conditions for the cusp are then given by d a Z = 0 
and d a X = 0 as discussed in Section |4.5.1|. The latter condition is automatically 


satisfied near the boundary because of d a X ~ a, while the former gives t 0 = 0. 
In Eq. (|4.27|) , the denominator has t 0 , and it appears that the numerator does not 
cancel the zero in the denominator. Hence, it is natural that the force diverges when 
the cusp arrives at the boundary. 


5.4 Z 2 -symmetric quench 

When the two endpoints of the string are simultaneously quenched in the opposite 
directions with the same amplitude, there is a new contribution from the one-sided 
quench that the propagating waves collide at the ^-symmetric point X = 0, and 
cusps are expected to form on the collision. This case is equivalent to impose the 
Neumann boundary condition at X = 0. Such a condition is typically imposed in 
probe D-brane embeddings. We would like to emphasize that understanding the 
difference between this case and the case without the Neumann condition would be 
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(a) 0 < t/L < 13 (b) 0 <t/L< 30 


Figure 13. Time dependence of the forces acting on the quark and the antiquark for the 
longitudinal one-sided quench with e = 1.0 x 10 -2 and A t/L = 2. For these parameters, 
cusps form on the string. Pulse-like oscillations in the forces are getting amplified as time 
increases and diverge eventually. 


important for distinguishing mechanisms for cusp formation whether the cusps are 
formed spontaneously as discussed in previous sections or formed with a help of the 
collisions or the Neumann condition. 

We investigate the cusp formation by computing the condition (|4.20| ). Results 
are shown in Fig. |h| for varying e with A t/L = 2 fixed. The left panel shows the 
times for ( [f.20|) to be satisfied for the first time, and the corresponding ^-coordinates 
are plotted in the right panel. We find that there are two kinds of cusp formation 
and wave collisions at the Z 2 -symmetric point: One is cusp formation by the wave 
collision, and the first cusp formation in this case is marked with red points in the 
plots, since in this case these cusps disappear once the colliding waves pass. The 
other is that cusps are formed on the propagating waves in the same way as in 
the one-sided quench, and the cusp formation for this case is marked with green 
triangles. The times for the cusp formation are clearly discretized, and the locations 
are concentrated to 0 = 0o/2. The small change in the formation time is because 
the propagating speed of the wave slightly varies for different e. 

Closely looking at the results in evaluating (|4.20|) , we find that even numbers of 
cusps are created for the first case (red points). In the bulk coordinates, the cusp 
formation is seen at very close to 0 = 0 o /2 and one might naively think that only one 
cusp was generated at a point. In the worldvolume results, however, when the waves 
collide at 0 = 0o/2, we find that a pair of cusps whose orientations are opposite 
are created, and then these cusps pair-annihilate shortly when the waves pass by. 
Hence, the cusps created by the collision do not propagate away from that point, and 
therefore these instantaneous cusps are not observed in the force at the boundary. In 
this case, other formation of cusps in the same way as that in the one-sided quench 
also happens afterward. 
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(a) Cusp formation times 



(b) Corresponding formation points 


Figure 14. Cusp formation times for the ^-symmetric quench with At = 2. Cusp creation 
by wave collisions around 0 = /3q/2 are plotted with red dots, while events apart from that 
point and in the same way as the one-sided quench are plotted with green triangles. 


In the cases of the green triangles in Fig. [L4|, the cusps are formed slightly before 
the collision point, and these cusps subsequently collide at the Z 2 -symmetric point. 
In fact, it is seen in Fig. |14(a)| that these points go ahead of the cusp formation by 


the collisions. These cusps then continue traveling on the string, inferring that the 
waves are already magnified enough for forming cusps. 

In the ^-symmetric case, it is convenient to evaluate the time evolution of the 
worldsheet Ricci scalar at 0 = /3o/2 since many events of the cusp formation occur 
there. The Ricci scalar is given by 


R = 


2 ip/uv,u r Juv,v r jux.Auv,uv) 


'■y3 
I uv 


(5.1) 


For the static configuration, this monotonically changes from R = — 2/T 2 at 0 = 0 
to R — —4/£ 2 at 0 = 0o/2, and because of this coordinate dependence it might 
be desirable to compare the Ricci scalar at a fixed 0. The Ricci scalar diverges 
on top of a cusp, and the cusp formation condition (|4.20|) is consistent with the 
divergence of the Ricci scalar since j uv in the denominator becomes zero when ( |4.20|) 
is satisfied. Practically, as we use discretized computations, the Ricci scalar does not 
exactly become infinity, while at least it becomes huge. Results of the Ricci scalar 
representing the cusp formation at second, third, and forth collisions are shown in 
Fig. [15]. In these results, the absolute value of the Ricci scalar suddenly becomes 


huge at the cusp formation, and come back to of order one when the waves pass. In 
our setup, the numerical evolution does not breakdown after the Ricci scalar diverges 
in contrast to the D3/D7 case [T0|—|T2| , although finite-iV c corrections may have to be 
taken into account. 
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Figure 15. Ricci scalar evaluated at <f> = /3 q/2 when A t/L = 2. 
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Figure 16. Snapshots of the string in the transverse linear quench. We set the parameters 
as A t/L = 2, e = 0.03. 

6 Results for the transverse linear quench 

6.1 Cusp formation 

In this section, we study the string dynamics induced by the transverse linear 
quench ( |4.4| ). With such a quench, the string moves in the (3+l)-dimensions spanned 
by (t, z, xi, x^). Snapshots of string configurations under a quench with parameters 
At/L = 2 and e = 0.03 are shown in Fig. |16[ The left panel is just after the quench, 
t/L = 2, 2.2, 2.4, where we do not find cusps. However, in the right panel for late 
time t/L = 15, 15.2, 15.4, we find cusps on the string. This demonstrates that cusps 
can form in the transverse linear quench. For quenches with smaller amplitudes 
(e < 0.02), we did not find cusp formation for the period we computed the time 
evolution. 

6.2 Energy spectrum in the non-linear theory 

We expect to see nonlinear origin for the cusp formation in the energy spectrum also 
in the case of the transverse linear quench. I 11 Fig. [L7|, we show the time dependence 
of the energy spectrum when the parameters are A t/L = 2 and e = 0.03. The 
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Figure 17. Time dependence of the energy spectrum in the transverse linear quench with 
A t/L = 2 and e = 0.03. Magenta points correspond to the time for the cusp formation. 
We fit them by e oc n~ a and show the result by a dotted line. 


time for the cusp formation is t/L = 14.45 for these parameters. Red points are 
just after the quench, and magenta points correspond to the time at cusp formation. 
We find the direct energy cascade as in Section pT2| . and eventually the spectrum 
obeys a power law until the time of cusp formation. Thus, also in the transverse 
linear quench, we find the turbulent behavior toward the cusp formation. Fitting 
the spectrum by £ oc n ~ a , we obtain a = —1.533 ± 0.133. The exponent has similar 
value as that for the longitudinal quench. 


6.3 Time-dependence of forces acting on the heavy quarks 


We turn to the time-dependence of the forces acting on the quark and the antiquark 
in the transverse linear quench. Since the motion of the string is in the (z,xi,X 2 )~ 
space, the X\- and x 2 -components of the forces can be non-zero. In Fig. |18(a)| , we 
show ( F) and ( F) as functions of time for A t/L = 2 and e = 0.01, with which cusps 
do not appear on the string. Similar to the longitudinal quench, pulse-like oscillations 
are repeated at intervals. Although there are sharp peaks, they are always 0( 1) in 
units of A -1/,2 L 2 . We also find that (Fi(f)) < 0 and (Fj(f)) > 0. Thus, the force 
between quarks is always attractive. 

In Fig. |18(b)| , we show the absolute values of the forces for At/L = 2 and 
e = 0.03, where cusps are formed as seen in Section |6T|. The pulses are getting sharp 
and amplified as time increases, and eventually after the cusp formation, the forces 
diverge when the cusps arrive at the boundaries. We also monitored to, which is 
in the denominator in Eq. ( |4.27| ), at the boundaries, and found that it is consistent 
with zero at t/L ~ 13.6 and 15.9 at x = x q and Xq, respectively. 
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Figure 18. Forces acting on the quark and the antiquark in the transverse linear quench. 
We fix the time scale of the quench as A t/L = 2. Left and right panels are for e = 0.01 
and e = 0.03, respectively. In the right panel, we take the absolute values of the forces, 
and the vertical axis is log scale. 


7 Results for the transverse circular quench 


Finally, we study the string dynamics induced by the transverse circular quench ( |4.5|) . 
The string moves in all (4+l)-dimensions spanned by (t, z, Xi,x 2 , £ 3 ). For the trans¬ 
verse circular quench, we did not find any cusp formation at least for modest pa¬ 
rameters: around at e ~ 0.01 and A t/L ~ 1. Nevertheless, we found an interesting 
behavior in the energy spectrum. In Fig. [19], we show the time dependence of the 
energy spectrum for parameters A t/L = 2 and e = 0.02. In the early time evolution 
until t/L < 14, there is a direct energy cascade: the energy is transferred from large 
to small scales, and eventually the spectrum obeys a power law at t/L ~ 14. Fitting 
the numerical data at t/L = 14, we obtain e n oc n - 2 019 ± 0 029 . For t/L > 14, however, 
we find that this turns into an inverse energy cascade: The energy is transferred to 
the large scale. Thus, in contrast to the previous low dimensional cases, the power 
law once realized at an intermediate time t/L = 14 is not maintained in the late 
time. 

In Fig. H we show snapshots of string configurations around the “turning point” 
of the energy cascade: t/L = 14,14.2,14.4. Since the string motion is in the (4 + 
l)-dimensions, we project the string profile into {xi,x 2 ,z)- and (aq, x 3 , z)-spaces. 
Although cusp-like points can be seen in the right figure, these are not real cusps: 
We find that although roots of Jz, Jx 2 and Jx 3 become close at t/L ~ 14, Jx 1 is not 
zero at the point. (See Eq. ( |4.20|) .) However, the perturbation variable y 3 defined in 
Eq. (|4.23| ) becomes cuspy, namely, its energy is transferred to the small scale. Hence 
the direct energy cascade appears in t/L < 14. After t/L ~ 14, the cuspy shape gets 
loose because of the dispersive spectrum in the linear perturbation. 

In Fig. EH], we show the time dependence of the forces acting on the quark and 
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(a) Direct energy cascade (b) Inverse energy cascade 


Figure 19. Time dependence of the energy spectrum for the transverse circular quench 
with A t/L = 2 and e = 0.02. 



Figure 20. Snapshots of the string in the transverse circular quench with A t/L = 2 and 
e = 0.02. In the left and right figures, we project the profiles of the string into (x\,X2,z)- 
and (xi, xs, z)-spaces, respectively. 

the antiquark for e = 0.02 and A t/L = 2. We do not find the divergence of the 
forces. However, the forces are magnified until t/L ~ 14 and can be repulsive at 
some time intervals ((ifi) > 0 and (iq) < 0) even though there is no cusp formation. 
The magnitude of the forces in the late time is not as big as that period, reflecting 
the looseness in the cuspy shape. 

8 Summary and Discussion 

We studied nonlinear dynamics of the flux tube between an external quark-antiquark 
pair in Af = 4 SYM theory using the AdS/CFT duality. We numerically computed 
the time evolution of the string in AdS dual to the flux tube when we perturbed the 
positions of the string endpoints to induce string motions. We considered four kinds 
of quenches that were chosen to represent typical string motions: (i) longitudinal one¬ 
sided quench, (ii) longitudinal Z 2 -symmetric quench, (iii) transverse linear quench, 
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t/L 


Figure 21. Time dependence of the forces acting on the quark and the antiquark for the 
transverse circular quench with e = 0.02 and A t/L = 2. 

and (iv) transverse circular quench. (See Eqs. (|1.1| - |4.5|) and Fig. |B|.) For (i)-(iii), we 
found cusp formation on the string. In the time evolution of the energy spectrum, 
we observed the weak turbulence, that is, the energy was transferred to the small 
scale, and the energy spectrum eventually obeyed a power law until the time of the 
cusp formation. The cusp formation occurred only when the amplitude of the quench 
was larger than a critical value, e > e cr i t , and the dependence of its magnitude on 
the quench duration At was given by a simple form e cr ; t oc (A t/L) 3 in small A t/L. 
When the cusps arrived at the AdS boundary, we observed the divergence of the 
force between the quark pair. For (iv), we found no cusp formation. Nevertheless, we 
observed a direct energy cascade and the power law spectrum for a while. However, in 
late time the direct cascade turned into an inverse energy cascade, where the energy 
was transferred to the large scale. There was no divergence of the force between the 
quark pair. 

How can we understand the weak turbulence of the string in view of gauge 
theory? Eigen normal modes e n of the fundamental string studied in section |3] can 
be regarded as the excited states | n) of the flux tube in the gauge theory side. Hence, 
the fluctuating string solution, such as R(t,(j>) = zo + L2 n c n(t)e n ((f)), corresponds to 

OO 

IV') = |0) + ^c n (f)|n) (8.1) 

n= 1 

in the boundary theory, where |0) is the ground state. The weak turbulence implies 
that |c n (t)| with n 1 tends to increase as a function of time. Therefore, in the 
late time, the probability of observing highly excited states is high compared to the 
linear theory. Similar phenomenon has been found in the D3/D7 system dual to 
J\f = 2 supersymmetric QCD [|T0|— , where a direct energy cascade was found in 
the fluctuations on the D7-brane and regarded as production of many heavy mesons 
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in the SQCD. In that paper, this phenomenon was referred as “turbulent meson 
condensation”. Although the endpoints of the flux tube we considered are regarded as 
nondynamical and infinitely heavy quarks, the string turbulence found in this paper 
would be regarded as the microscopic picture of the turbulent meson condensation. 

We found cusp formation when the motion of the string is restricted in (2 + 
1)- and (3 + l)-dimensions. The divergence of the forces acting on the quarks is 
accompanied by the cusp singularities, where we expect that finite-1V C effects will 
become important. These will contain quantum effects of the string, and such effects 
may resolve the cusp singularities and the divergence of the forces. Nevertheless, the 
cusp formation in the classical sense can give us observable effects: Finite- N c effects 
will also appear as gravitational backreactions. If these are taken into account, a 
strong gravitational wave will be emitted at the onset of the cusp formation!] For 
cosmic strings in flat spacetime, gravitational wave bursts from cusps have been 
studied in Ref. |42|], and it has been found that their spectra obey a power law in the 
high-frequency regime. It would be nice to compute the gravitational waves from the 
string in AdSs and find the power law spectrum. The description in the dual field 
theory may be the power law spectrum in gluon jets from the flux tube. 

When the motion of the string is in (4 + l)-dimensions, we did not find cusp 
formation. Hence, in AdSs spacetime, the cusp formation on the string is not a 
general phenomenon but accidental one. This implies that the dual phenomenon 
to the cusp formation is not ubiquitous in the (3 + l)-dimensional boundary field 
theory. However, if we consider the many-body system of quark-antiquark pairs, it 
is possible that the time evolution of some flux tubes happens to be restricted in 
lower dimensional spaces approximately, Q and such flux tubes would be able to emit 
the gluon jets with the approximately power law spectrum, which is characteristic 
to the cusp formation. Besides, gravitational wave bursts in the presence of extra 
dimensions were discussed in [[43], fh|]. Even though real cusps are not formed, there 
may be some gluon emission from cuspy shapes. 

There are some future directions in our work. In this paper, we only considered 
modest values for the amplitude of the quench, e ~ 0.01. For a large value of e, 
we expect that the string can even plunge into the Poincare horizon because of the 
strong perturbation. This will demonstrate a non-equilibrium process of breaking 
of the flux tube. It is also straightforward to take into account Enite temperature 
effects in this process. It will be also interesting to consider the string motions 
in confined geometries [jT7| . In the theories dual to these backgrounds, the quark- 
antiquark potential is linear, and the presence of such potential may affect conditions 
for cusp formation. Studying nonlinear string dynamics in such geometries may give 


9 In the asymptotically flat spacetime, gravitational self-interaction of cosmic strings has been 
perturbatively studied in Ref. pl|. This work suggests that cusps survive the backreaction. 
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We thank Claude Warnick for pointing out this argument. 






new insights into understanding the QCD flux tubes and non-equilibrium processes 
in realistic QCD. 


Closed strings rotating in AdS and having cusps were constructed in Ref. [05 


Although our dynamical cusp formation on an open string is different from the ex¬ 
istence of cusps in those steady solutions, it may be interesting to obtain useful in¬ 
formation from such configurations. In JT0|, the universal exponent in the power law 
was deduced from a stationary solution called critical embedding in the D3/D7-brane 
system in the presence of a constant electric field, and results in time dependent com¬ 
putations obeyed that universal value. In our setup, we do not have a corresponding 
static cuspy configuration, but our power law exponents, distributed around 1.4, may 
be naturally understood from cuspy stationary strings in AdS. 

Ultimately, it will be important to understand the mechanism relevant for the 
the turbulent behavior. For the integrable Wilson loops such as those in AdS 5 X S 5 , 
the turbulent behavior may be studied with the techniques of integrability. 

In non-linear systems, there may be underlying chaos. In [[fi|, closed strings 
moving in Schwarzschild-AdS background were studied from the viewpoint of chaos. 
It may be interesting if chaos is seen also in the motion of open strings and the 
turbulent behavior is understood, particularly in non-integrable situations. 
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A Numerical methods 


In this appendix, we explain our numerical method for solving the equations of 
motion of the string ( |4.11| ). Basic ideas are explained in Appendix A in 0. 

We found that the original form of the equations of motion ( |4.11 ) is numeri¬ 
cally unstable. To stabilize the numerical evolution, it is effective to use the con¬ 


straints (0.81) . From them, we have 


T 


, u 


= jzi + x%, 


T 


,V 


= JZ% + X% , 


(A.l) 
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Figure 22. Discretization of the world volume. 


where we choose the positive signs for the square roots since we take d u and d v as 
future directed vectors. Eliminating T u and T jV from Eq. ( |4.11|) , we obtain 


T„v = j[(Z\ + X%)WZ,„ + (Z 2 „ + X%)'VZ ,„] , 

Z,„„ = i[(Z, 2 „ + Xlfl\Zl + X 2 ,,) 1 ' 2 + Z,„Z,„ - X„ ■ X,,] , (A.2) 

X tUV — ~^(X :U Z iV + X >V Z >U ) . 

The evolution equations in these expressions are found numerically stable. 


To numerically solve (|A.2|) , we discretize the world volume (w, w)-coordinates 
with the grid spacing h as shown in Fig. [22|. Let us denote the fields (T, Z , X) by T. 
At a point C apart from the boundary, the fields and their derivatives are discretized 
with second-order accuracy as 






N 


— T 


w 




uv\C 


Me = 




N 


h? 

V Lr — $«/ 


Me = 




N 


Tr + T 


w 


4b 


2/i 




Tr + T 


T| c = . (A.3) 

2 h 1 |G 2 v ; 

Discretization error is 0(h 2 ). Substituting these into the evolution equations (|A.2|) , 
we obtain nonlinear equations to determine Tat by using known data of Tr, Tw, 
and Tg. We use the Newton-Raphson method for solving the coupled nonlinear 
equations. 

The equations for the boundary time evolution ( |4.15|) become coupled nonlinear 
equations of T,y and X^. ( Z^ = 0 is trivially imposed.) The form at a = 0 is 


Tn = Ts 


4 Z 2 W + {X N -X 


Xv = x q (T 1 


N 


(A.4) 


where we used a relation for Z, Z E — — z w , derived from the boundary condition 
Z tUV = 0. At the other boundary a = /3q, Zw and x q in (|A.4|) are replaced with Z E 
and Xq. These equations are also solved by using the Newton-Raphson method. 
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Figure 23. Constraint violation for several resolutions, N = 200, 400, 800, 1600. (a) Lon¬ 
gitudinal one-sided quench with A t/L = 2 and e = 0.01. (b) Transverse linear quench with 
A t/L = 2 and e = 0.03. (c) Transverse circular quench with A t/L = 2 and e = 0.02. For 
(a) and (b), the right ends of the figures correspond to the times for the cusp formation. 


B Error analysis 

In this section, we estimate errors in our numerical calculations. We define 

Cl = -0-r 2 „ + z% + x%) , c 2 = T(_rj + z 2 v + X%) . (b.i) 

These constraints should be zero for exact solutions. Hence, these can be nice indi¬ 
cators of our numerical errors. For visibility of the constraint violation, we introduce 

CmaxO) = max(|Ci|, \C 2 \) , (B.2) 

fixed v 

where we take the maximum value when we vary u on a fixed v surface. We also 
choose the bigger of the two constraints, \C\\ and | C 2 \. Introducing an integer N 
such that the mesh size is given by h = (3 0 /N, we plot C max (u) for several values 
of N in Fig. [23|. We see that the constraint violation is small ( C mSLX ~ 10~ 3 even 
for N = 200) and behaves as C max oc 1/N 2 . This is consistent with the fact that 
our numerical method has the second order accuracy. In this paper, we mainly set 
N = 800. Then, the constraint violation is 0(1O -4 ). 


C Energy spectrum in the linear theory 

In this appendix, we derive the energy spectrum induced in the linear theory of 
section when a quench is added on the boundary. The equations of motion for 
the perturbation variables (yi, y 2) X 3 ) are given in Eq. ( |3.8|) . Here, we focus on the 
longitudinal mode xi f° r simplicity and denote it as xi — X■ Application to trans¬ 
verse modes is straightforward. For the quench, we consider the following boundary 
conditions for y: 

x(t, 0 = 0)= Xb(t) , x(t, 0 = A)) = 0 , (C.l) 
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where Xbif) is the quench function assumed to have a compact support at 0 < t < At. 
We also assume that the solution is trivial before the quench, x(t < 0, 0) = 0. 

We firstly consider a time independent solution to ( |3.8|) , 


US{(j>) = 0 . 


(C.2) 


A solution is given by 


where A = f^° dcj)'/ h(4>') is a constant for normalization, 
function S behaves as 


s(d>) = - 


A 


r+ dcf)' 

4, KV) ’ 


5 = 1 — 


3Ar§ 


+ ••• (0~O), 5 = 


3Ar§ 


(A) - (Pf 


Let us introduce the quench Xb(f)- Using the function 

x{t, 4>) = x(t, 4>) + Xb(t)S((p). 


(C.3) 

Near the boundaries the 

+ ••• 0£~A)). (C.4) 

S, we define x as 

(C.5) 


The new variable x satisfies trivial boundary conditions, x(t, 0 = 0)= y(t, 0 = 0 O ) = 
0. The equation of motion for x becomes 


{d 2 t + U)x = —Xb(t)S(<f>) , (C.6) 


where we used Eq. (|C.2| ). 

To solve the equation, we consider a Green’s equation: 

{d 2 t + n)G(t, f- 0 , 0 ') = s(t - H)s (0 - 0 '), 


(C.7) 


where G is the Green’s function. By using G, a special solution to Eq. ( |C.6|) can be 
written in the form 

/ oo 

dt' / d4)'G(tA'-A,<j>')Xb(t')S{4>') ■ (c. 8 ) 

-oo J 0 

Operating dt on both sides of the Green’s equation and taking the limit of 
e —y 0, we obtain the junction condition as 


d t G | 


t=t '+o 
t=t'—0 


<5(0 - 00 . 


(C.9) 


For t < t', we assume the Green’s function is trivial: G — 0. For t > t', the Green’s 
function can be written as 


G = ^2[a n smuj n (t - t’) + b n cosu n (t - t')]e n {4>) , 

n 


(C.10) 
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From the continuity of G at t — t' we have b n = 0, and then from the junction 
condition (|C.9|) , we find that a n satisfies 


^ ^ d n ,cu n e n (0) b(cp cj )) 


(C.ll) 


Operating (e n , *) to the above equation, we obtain 


a n = —7 (</>')e n (<//) 
(jJ n 

Thus, the Green’s function can be written as 


(C.12) 




o 


(t < t '), 

En W n lsinW n(i-^)7(0O e n(0') e n(</>) 0 > t') . 


(C.13) 


Since the Green’s function is zero at t < t', the special solution obtained from 
Eq. (|G. 8 |) is also zero before the quench, t < 0, and this is nothing but the solution 
we are looking for. After the quench t > T, the solution becomes 


fT 

X = - y2u>- 1 S n e n (<f>) / dt'xo{t') Sin u n (t - t') 
r T 

= ^2uj n S n e n ((p) / dt'xb(t')smuj n (t-t') , 


(C.14) 


where S n = (S,e n ). Note that x = X after the quench. At the second equality, we 
integrated by parts twice. 

It is then straightforward to compute the energy spectrum. The mode coefficient 
c n = (x, e n ) is computed by using (|C.14|) as 


c n (t) = uj n S n / dt’xbit’) sin u n (t - t ') , 


(C.15) 


and then from Eq. ( f4.25|) the energy spectrum for the longitudinal quench in the 
linear theory becomes 


\/A^o r- 


47T 


K + W « C n] = 


VXzouiSl 


47T 


xM\ 2 ■ 


(C.16) 


where we define x(u>) = f dtXb(t)e~ lut . The energy spectrum does not depend on 
t as we expect. Taking into account the transverse modes, we obtain the energy 
spectrum as 


£n. 


VXZO 

An 


UJ n S l\xM\ 2 + uj' n 4 S' n 2 \xiM\ : 

i=2,3 


(C.17) 
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where 


XiM = / dtXb(t, <fi — 0)e 


—iujt 


s' = + r ° , 5;=(S',e„)', (C.18) 

Zl n 


and Xfe are the quench functions for Xi- The total energy is e = Y^=i £ n- 

The spectrum (|C.17|) is dehned only for integer n. However, in Figs. 00 and 
P~9| , we generalize e n to a continuous number by interpolating co n , S n , u>' n , and S' n and 
regarding them as function of continuous number n for visibility. 


D Forces acting on the quark and the antiquark 

In this appendix, we derive the formula for the forces acting on the quark endpoints 


( 14.27|) . We denote the on-shell Nambu-Goto action as S[x q ,x q ], where x q and x q 
are the locations of the string endpoints regarded as the quark and the antiquark, 
respectively, at the AdS boundary: X(t, z —> 0) = x q , x q . The on-shell action relates 
to the partition function of the boundary theory as 


Z C FT[x q , x q ) = e iSlx *’ x * ] . 
In the field theory, the partition function is written as 


(D.l) 


Z C ft [x q , Xg\ = J V(j) exp (iS'sYMM + iS q y>(x q ), x q ] + iS q [(j)(x q ), x q ]) , (D.2) 

where 0 represents the set of the fields in Af = 4 super Yang-Mills theory and 
Ssym is its action. S q and S q denote the actions for the quark and the antiquark, 
respectively. We regard x q (t) and x q {t) as external fields. The quark action is 
schematically written as 


^q\(fo(xg] : Xg\ / dt 


-m 


y/l — v 2 + L 


int yPyXq) i X q 


(D.3) 


where m is the quark mass and L- mt corresponds to the interaction term with SYM 
fields, and the velocity of the quark is introduced as v = dx q /dt. We define the force 
acting on the quark as 

F = Sx~ J dtLint ^ x i^ x <^ ■ ( D - 4 ) 

The variation of the on-shell Nambu-Goto action is then related to the field theory 
terms as 

A ,<? A 

(D.5) 


ss s , „ . , 

CFT = ~ mhvy + { } ’ 


where (F) = Z~ x f V(f)Fe l ‘ s s YM + i S q +iSq ^ ^ _ \jy/\ _ and the dot in the upper 
right of the parentheses denotes ' = d/dt. In the first equality, we used the AdS/CFT 
duality ( p.l|) . 
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Now, we evaluate 5S/5x q in the gravity side. For this purpose, it is convenient to 
use target space coordinates (t, z ) as the world sheet coordinates. Then, the position 
of the string is specified by x = X(t, z ), and the Nambu-Goto action becomes 


S = -0 J dtdz 4 [(1 - ^ 2 )(1 + X' 2 ) + (X • X') 2 ] 1/2 . (D.6) 

The equations of motion are given by 


’-(1 + X ,2 )X + (X • X')X'" 

1 

(1-X 2 )X'+ (X-X')X 


1 

Z 2 Vl 


where = d/dt, ' = d/dz, and £ = (1 — X 2 )(l + X' 2 ) + (X ■ X') 2 . Solving these 
near the AdS boundary, we obtain 

X = x q (t) - \p 2 az 2 + x 3 ^ 3 + @( z4 ) ) (D.8) 

where a = d?x q /dt 2 . Let us consider the variation of the on-shell action with respect 
to one of the endpoints of the string, x q , and the other endpoint is fixed, 8x q = 0. 
By defining Lagrangian as S = f drdaC , the variation of the action becomes 


8S[x q ,x q ] = / dtdz 


S(d a X) 


dC 


= i dtdz a a [sx- g{ ^ x 


d(d a X 

dC 


8X ■ 


dC 

ax 


dt dX ■ 


dC 


dX' 


(D.9) 


X=Xq,Z=e 


where we take the cutoff at z = e. At the second equality, we used Euler-Lagrange 
equation, dC/dX = d a [dC/d(d a X)\. There is no contribution from the other bound¬ 
ary X = x q and z = e since 5X = 0 there. Substituting Eq. (p.8|) into Eq. (p.9|) , 
we obtain 


= J dtSx,- + ^~( x ^ + 7 2 (»’' ® 3 )v) 


(D-10) 


Hence, the upshot for dS/dx q is 
dS 


dx r , 


2ne 


(7 V Y + 7T-^( X 3 + 1 2 ( v ■ x 3)v) . 

Z7T7 


(D.ll) 


Comparing above expression with Eq. (P.5|) , we can see that the hrst term in ( D.ll| ) 
corresponds to a diverging quark mass m ~ 1/e. This is a natural consequence since 
we are considering an infinitely extended string. Setting m = -\/A/(27re), we obtain 
the force acting on the quark from Eqs. ( p.5|) and (P. 11|) as 


i F i(t)) = 0( S v +7 2 ^i)^X j | z=0 , 


(D.12) 
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where we replaced a? 3 with 6^X\ z= q/§. Note that there can be a finite difference 
between m and \/A/(27re), and an extra-term proportional to ( 71 ;)' may appear in 
Eq. (p,12|) . However, it can be eliminated by adding a local counter term propor¬ 
tional to \/l — v 2 in the quark action S p . The same formula can be applied to the 
force acting on the antiquark at Xq. 
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